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Abstract 

The  Aero  Propulsion  Laboratory  is  currently  developing 
a  computer-aided  design  program  for  high  power  airborne 
systems.  An  important  part  of  this  design  program  will  be 
the  feasibility  study  which  was  to  be  based  on  summary  algor¬ 
ithms.  These  algorithms  were  to  relate  the  weight  and  volume 
of  each  system  component  to  the  contributing  operating  para¬ 
meters  . 

This  study  first  centers  on  the  requirements  for  an 
interpolation  scheme  to  form  the  summary  algorithms.  The  scheme 
will  need  to  work  in  at  least  four  dimensions  and  produce  accu¬ 
rate  results  over  wide  ranges.  The  advantages  of  the  two 
possible  interpolation  approaches  (algorithm  development  and 
direct  interpolation)  are  described.  The  results  of  this  com¬ 
parison  show  clear  advantages  in  the  direct  interpolation 
approach  using  stored  data. 

The  remainder  of  this  study  is  a  description  and  evalu¬ 
ation  of  the  subroutine  which  was  developed,  INTERP.  It  has 
the  flexibility  to  interpolate  a  data  array  with  two  or  more 
independent  variables  and  vTill  output  the  values  of  any  number 
of  dependent  variables.  The  routine  also  compensates  for  mis¬ 
sing  values  in  the  known  data  array  and  issues  an  error  message 
to  the  user  when  a  test  point  beyond  the  data  range  is  input. 
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INTERPOLATING  SUBROUTINE  FOR 


HIGH  POWER  SYSTEM  DESIGN 

I.  Introduction 


Background 

The  Aero  Propulsion  Laboratory  has  responsibility  for 
the  development  of  multi-megawatt,  airborne  power  systems. 
These  high-power  systems  will  be  required  for  directed  energy 
weapons  being  developed  at  the  Air  Force  Weapons  Laboratory. 
The  components  required  in  these  power  systems  will  include 
turbine,  generator,  transformer,  rectifier,  and  pulse¬ 
forming  network.  The  three  steps  in  the  system  development 
are  (1)  the  system  feasibility  study,  (2)  the  detailed  com¬ 
ponent  design,  and  (3)  the  dynamic  system  simulation 
(Figure  1) . 

The  first  step,  system  feasibility  study,  will  use 
summary  algorithms  which  relate  the  weight  and  volume  of  each 
component  to  the  contributing  operating  parameters.  These 
algorithms  will  then  be  applied  over  broad  parameter  ranges 
to  determine  optimal  combinations  of  components.  Operating 
parameters  determined  in  this  step  will  then  be  converted  to 
component  design  specifications.  Computer  programs  which  can 
produce  detailed  component  designs  will  then  be  used  in  step 
two  to  specify  the  detailed  design  for  each  component. 

The  final  step  will  be  the  dynamic  computer  simulation 
of  the  inter-connected  components.  This  simulation  will 
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Fig  1.  Flowchart  of  Computer  Aided  Design  Concept 
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evaluate  electrical  and  thermal  performance  of  the  entire 
system.  The  electrical  simulation  will  permit  evaluation  of 
voltage  and  current  transients  and  aid  in  the  refinement  of 
the  control  system.  Thermal  simulation  will  provide  the  data 
needed  to  specify  the  cooling  system  requirements. 

The  system  study  is  still  in  the  first  phase,  the 
development  of  summary  algorithms.  A  summary  algorithm  des¬ 
cribes  the  weight  or  volume  of  a  component  as  a  function  of 
operating  parameters.  Examples  of  operating  parameters  are 
input  voltage,  frequency,  and  output  power. 

Currently,  the  Aero  Propulsion  Laboratory  is  using 
civilian  contractors  to  develop  design  programs  and  point 
designs,  which  cover  the  expected  range  of  each  of  the  para¬ 
meters.  One  contract  is  with  the  Raytheon  Corporation  to 
provide  120  minimum  weight  and  volume  designs  for  rectifiers. 

The  contract  also  includes  110  inverter  designs  and  65  designs 
for  pulse- forming  networks.  These  designs  will  be  delivered 
during  1981  and  1982. 

Problem 

The  current  problem  is  the  need  for  a  computer  pro¬ 
gram  which  can  use  the  data  provided  in  these  designs  to  form 
the  summary  algorithms  needed  for  system  feasibility  analyses. 

The  first  objective  of  this  study  was  to  determine  the 
best  format  for  these  algorithms.  The  format  will  be  deter¬ 
mined  by  the  choice  of  a  multi-dimensional  interpolation  scheme. 
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Considerations  in  this  choice  will  include  (1)  compati¬ 
bility  of  the  algorithm  format  with  the  next  step  of  the  system 
analysis,  (2)  compatibility  of  the  scheme  implementation,  and 
(3)  complexity,  range,  and  accuracy  of  the  scheme.  The  actual 
writing  of  the  computer  program  to  form  these  algorithms  will 
be  the  second  part  of  this  study. 

Approach 

The  first  step  in  this  investigation  was  a  review  of 
the  High  Power  Studies  (a  set  of  technical  reports  prepared 
for  the  Aero  Propulsion  Laboratory)  and  direct  consultation 
with  project  engineers  in  the  laboratory  to  consider  the  ques¬ 
tion  of  compatibility  with  the  analysis  program.  The  other 
area  of  investigation  was  a  literature  search  and  consultation 
with  Dr.  Lee,  Chairman  of  the  AFIT  Mathematics  Department,  to 
compare  interpolation  schemes. 

Sequence  of  Presentation 

A  more  detailed  analysis  of  the  requirements  which 
affect  the  choice  of  the  interpolation  scheme  is  presented  in 
Chapter  II.  Chapter  III  describes  the  two  possible  approaches 
to  the  interpolation  problem  and  the  advantages  of  each.  The 
development  of  an  interpolating  subroutine  based  on  the  direct 
interpolation  scheme  is  also  proposed. 

Chapter  IV  describes  the  subroutine  developed,  INTERP. 

In  Chapter  V  the  accuracy  and  effectiveness  of  INTERP  is  evalu¬ 
ated  using  the  program  CGEN  and  stored  data  based  on  conventional 


generator  designs.  Conclusions  and  recommendations  for  future 
study  are  presented  in  the  final  chapter. 

There  are  also  two  appendices.  Appendix  A  includes  a 
"User  Guide"  to  INTERP  along  with  a  sample  program  and  a 
listing  of  the  subroutine.  Appendix  B  contains  a  listing  of 
CGEN  as  well  as  some  of  the  inputs,  outputs,  data,  and  graphs 
used  in  the  evaluation  of  INTERP. 
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I I .  Investigation  of  Requirements 


As  stated  earlier,  the  primary  application  of  the  sum¬ 
mary  algorithms  will  be  in  the  computer  program  used  to  deter¬ 
mine  the  optimal  configuration  of  system  components.  This 
optimal  design  will  include  some,  but  not  necessarily  all,  of 
the  possible  system  components  indicated  below  with  their 
associated  parameters: 

1)  Turbine  (including  gas  generator  and  fuel  supply) -- 
speed,  power,  run-time,  number  of  starts 

2)  Generator  (three  different  types  will  be  considered) 
--speed,  voltage,  power 

3)  Transformer  --  voltage  in,  voltage  out,  power, 
frequency 

4)  Rectifiers  --  power  out,  frequency,  voltage  out 

5)  Inverters  --  voltage  in,  voltage  out,  power  out, 
frequency 

6)  Pulse- forming  networks  --  power,  pulse  width, 
pulse  repetition  time. 

Figure  2  depicts  a  possible  system  design. 

Summary  algorithms  -have  already  been  developed  for 
transformers,  using  careful  analysis  of  a  large  number  of 
smoothly  varying  design  points.  Some  algorithms  have  been 
developed  for  generators.  Due  to  discontinuities,  they  are 
valid  only  over  a  limited  range.  Below  is  an  example  of  one 
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Fig  2.  Generalized  High-Power  Airborne  System 

of  these  algorithms  (Ref  2:4). 

Specific  Weight  (Ibs/kW)  *  .157  [1.28  -  .28 
.  [-.06  +  •  [.8567  +  .1433(^)]  (1) 

where 

P  =  generator  output  power  (megawatts) 

RPM  =  rotational  speed 

V  =  generator  output  voltage  (volts,  line-neutral) 
Algorithms  developed  so  far  for  power  conditioning  components 
are  based  on  preliminary  designs  and  give  unacceptable  results 
in  most  cases. 

New  design  points,  now  being  developed  by  Raytheon, 
will  probably  provide  the  most  severe  test  of  the  interpolating 
scheme.  Large  data  sets  over  wide  ranges  will  be  developed. 
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Inverter  designs  based  on  four  independent  parameters,  with 
some  discontinuities,  will  be  included.  This  will  require  an 
interpolation  scheme,  in  at  least  four  dimensions,  which  can 
produce  reasonable  accuracy  over  wide  ranges. 
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III.  Theoretical  Approaches 

There  are  two  possible  approaches  to  an  interpolating 
problem.  The  first  is  to  use  the  available  data  to  generate 
an  algorithm,  similar  to  Eq  (1),  which  can  then  be  used  inde¬ 
pendently  of  the  data.  The  alternative  approach  is  to  retain 
the  data  in  an  array  or  file.  When  a  value  of  the  function 
is  needed,  it  is  necessary  to  find  the  nearest  known  points 
in  the  array  and  estimate  the  value  of  the  function  at  the 
new  point.  The  next  step  is  to  compare  these  two  approaches 
and  determine  which  is  best  suited  to  this  application. 

Algorithm  Approach 

First  to  be  considered  is  the  algorithm  approach.  This 
approach  has  been  used  to  date  by  engineers  working  on  the  high 
power  systems  under  consideration.  It  has  the  obvious  advan¬ 
tage  of  succinctness.  When  incorporated  into  the  overall 
design  feasibility  program,  the  algorithm  approach  allows 
faster  execution  and  requires  much  less  core  memory.  However, 
writing  a  computer  program  to  reduce  a  multi-dimensional  array 
of  data  to  a  compact  algorithm  proves  to  be  a  very  complex  task. 

One  way  to  design  such  a  program  would  be  to  use  the 
multi-dimensional  spline  technique  discussed  by  de  Boor  in 
Reference  3.  At  first  this  approach  appears  reasonable, 
though  quite  involved.  A  problem  appears  near  the  end  of  the 
development  where  a  system  of  perhaps  several  thousand 
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simultaneous  equations  must  be  solved.  Even  though  the  system 
has  considerable  structure,  so  that  highly  efficient  solution 
schemes  could  be  used,  this  is  a  difficult  problem.  An  alter¬ 
native  to  the  spline  approach  is  to  use  an  application  of  the 
multi-dimensional  Taylor  series  expansion  (Eq  2). 
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where 

(xj^,  X2,  ...  x^)  =  point  where  the  value  of  the 

function  is  known 

The  key  consideration  in  this  approach  is  the  "higher  order 
terms."  For  a  nonlinear  function  there  is  no  limit  to  the 
number  of  higher  order  terms,  each  of  which  becomes  more  com¬ 
plicated.  The  functions  being  considered  will  definitely  not 
all  be  linear.  However,  a  function  can  sometimes  be  made 
linear  using  a  change  in  variables.  For  instance,  if  a 
function  F  varies  as  x  +5,  we  could  substitute  X'  =  x  +5; 
and  F  would  be  linear  with  respect  to  x'.  This  is  the  approach 
which  has  been  used  on  several  data  sets  in  the  past.  For  data 
which  we  can  examine  and  manipulate  by  hand  with  the  help  of 
human  intuition,  we  have  a  reasonable  likelihood  of  discover¬ 
ing  variations  like  ax^  +  b,  ax^  +  b ,  e^  +  a,  or  a^  +  b. 


Physics  may  also  suggest  these  groups.  However,  to  program 
a  computer  to  uncover  all  of  the  possible  variations  would  be 


a  gigantic,  if  not  impossible,  task.  An  additional  problem 
is  that  these  algorithms  are  based  on  the  behavior  of  the 
function  near  the  midpoint  of  the  data  sets  and  often  produce 
serious  inaccuracies  when  applied  to  points  near  the  edges 
of  the  data  array. 

Direct  Interpolation 
Appro  aTE 

After  considering  the  complications  associated  with 
the  algorithm  approach,  a  careful  look  at  the  direct  inter¬ 
polation  approach  is  definitely  warranted.  An  advantage  to 
this  approach  is  that  output  is  limited  to  the  range  between 
the  neighboring  data  points.  This  avoids  the  possibly  extreme 
errors  of  a  modeling  approach.  If  more  accuracy  is  required, 
a  spline  or  other  more  powerful  interpolation  technique  can 
be  used. 

As  mentioned  previously,  this  approach  has  the  disad¬ 
vantage  of  requiring  more  time  and  memory  when  being  executed 
in  the  system  feasibility  program.  The  objective  of  the 
feasibility  program,  however,  is  not  to  run  in  minimum  time 
or  with  minimum  memory.  Whether  it  requires  ten  minutes  and 
15K  of  memory  or  six  hours  and  150K,  in  either  case,  if  the 
program  is  successful,  it  will  save  weeks  or  months  of  research 
and  greatly  reduce  the  acquisition  time  of  the  final  system. 

Another  disadvantage  to  be  considered  is  the  require¬ 


ment  to  store  and  maintain  the  data  sets.  This  is  normally 
done  on  a  tape  at  the  computer  center  and  requires  careful 
documentation  to  ensure  that  the  tapes  are  retained  and  that 


the  data  on  them  is  readily  accessible.  The  other  advantages 
to  be  gained  certainly  warrant  the  small  additional  effort 
required  to  document  these  data  sets. 

Conclusion 

An  objective  of  the  program  developed,  using  this 
interpolating  scheme,  will  be  its  applicability  to  operate  on 
data  sets  still  under  development.  To  meet  this  requirement 
and  to  ensure  reasonable  accuracy  from  the  scheme,  it  is  con¬ 
cluded  that  the  direct  interpolation  approach,  using  stored 
data,  is  the  best  for  this  application. 

The  next  step  in  the  study  is  the  development  of  a 
general  subroutine  to  accommodate  at  least  five  variables. 

This  subroutine  should  be  general  enough  to  be  used  on  a  data 
set  from  any  of  the  high  power  system  components.  Inputs  to 
this  routine  would  include  the  number  of  independent  vari¬ 
ables,  the  number  of  data  points  for  each  variable,  the  values 
of  the  variables  at  the  known  data  points,  an  array  containing 
the  function  values  at  the  data  points,  and  the  point  at  which 
the  function  is  to  be  evaluated.  An  advantage  to  this  scheme 
is  that  three  values  (weight,  volume,  and  efficiency)  can  be 
stored  in  the  same  data  array,  and  all  three  can  be  determined 
with  one  pass  through  the  subroutine. 
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IV.  Description  of  SUBROUTINE  INTERP 


SUBROUTINE  INTERP  was  developed  by  the  author  to  meet 
the  requirements  discussed  in  the  previous  chapter.  The  sub¬ 
routine  can  accommodate  two  or  more  independent  variables  and 
is  not  limited  to  five.  It  can  store  any  number  of  dependent 
values  at  each  point  in  the  array.  This  allows  storage  and 
interpolation  of  weight,  volume,  efficiency,  and  other  data. 

The  routine  includes  an  error  code  output  which  sig¬ 
nals  the  user  when  he  has  requested  data  from  a  point  outside 
the  stored  array  bounds.  If  data  are  unavailable  for  some 
points  in  the  array,  the  user  fills  these  locations  with  zeros 
and  the  routine  automatically  searches  until  non- zero  values 
are  found  on  which  to  base  the  interpolated  values.  Since 
these  values  normally  represent  values  of  weight,  volume,  or 
efficiency,  a  value  of  zero  should  not  otherwise  appear. 

In  the  following  section  of  this  chapter  is  a  dis¬ 
cussion  of  the  inputs  to  the  routine  and  how  the  data  arrays 
are  stored.  Subsequent  sections  give  an  overview  of  the 
routine  and  then  describe  each  of  the  three  main  sections  in 
more  detail.  A  complete  listing  of  INTERP,  a  user  guide, 
and  a  sample  program  are  included  in  Appendix  A. 

Arguments 

The  first  input  to  the  routine,  NDIM,  is  an  integer 
which  specifies  the  number  of  independent  variables  or 
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dimensions  to  be  handled.  Allowable  values  are  two  or  greater. 
The  second  input,  LV,  is  a  vector  (length  equal  to  NDIM}  which 
contains  integers  corresponding  to  the  number  of  data  points 
in  each  dimension.  These  points  do  not  need  to  be  evenly 
spaced.  The  order  of  these  points  establishes  the  order  in 
which  the  dimensions  will  be  considered  throughout  the  routine. 

LVT  is  the  next  input  and  is  the  sum  of  the  elements  in 
the  length  vector,  LV.  LDI  is  the  number  of  dependent  values 
stored  at  each  point  in  the  array.  Any  positive  integer  may 
be  used.  LDJ  is  the  product  of  the  elements  in  vector  LV. 

LVT,  LDI,  and  LDJ  are  used  along  with  NDIM  to  dimension  the 
remaining  arrays. 

The  next  input  is  the  vector  VAR  with  length  equal  to 
LVT.  The  values  of  the  independent  variables  at  the  known 
data  points  are  stored  in  this  vector.  The  values  of  the 
first  variable  are  stored  first  in  increasing  order,  followed 
immediately  by  the  values  of  the  next  variable  stored  the  same 
way,  and  likewise  through  the  independent  variables.  These 
variables  must  be  stored  in  the  same  order  as  the  length  values 
stored  in  LV.  The  input  WVE  is  the  array  containing  the  depen¬ 
dent  variables  (typically  weight,  volume,  and  efficiency). 

The  array  has  dimensions  LDI  by  LDJ.  The  first  index  of  the 
array  specifies  which  dependent  variable  to  consider  for  a 
given  set  of  independent  variables.  The  second  index  walks 
through  the  possible  data  locations  by  first  incrementing 
through  all  values  of  the  first  independent  variable,  then 
through  all  possible  combinations  of  the  first  and  second 
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independent  variables,  and  so  on  through  all  the  possible  com¬ 
binations  of  independent  variables.  Zeros  should  be  stored  at 
any  indices  where  the  data  are  not  available.  The  routine  will 
check  for  zero  values  and  search  for  the  next  non- zero  value. 

The  input  XV  is  a  vector  (length  equal  to  NDIM)  which 
contains  the  coordinates  of  the  test  point,  the  point  in  the 
data  array  where  the  values  of  the  dependent  variables  are  to 
be  estimated.  The  vector  LWK  (length  equal  to  NDIM)  is  a  work 
space  providing  temporary  storage  locations  for  index  values 
within  the  routine.  OUT  is  the  output  vector  of  length  LDI 
containing  the  computed  values  of  the  dependent  variables. 

The  final  argument  of  the  subroutine,  lERR,  is  a  two- 
digit  error  message  output.  The  least  significant  digit  will  be 
either  zero  or  nine.  Zero  indicates  that  a  requested  point  is 
less  than  the  lowest  point,  and  nine  indicates  a  requested  point 
greater  than  the  highest  data  points.  The  leading  digit  indi¬ 
cates  the  dimension  in  which  the  out  of  range  condition  occurred. 


Mathematical  Basis 

The  mathematical  basis  for  INTER?  is  the  multi¬ 
dimensional  Taylor  series  expansion  introduced  in  Chapter  III 
and  repeated  below  in  Eq  (2). 


F(x^,  X2,  ...  x^)  F(x^,  x^,  ...  x^) 

*  H,  (='1  ■  *  If,  (’'2  -  =<2) 


gp 

.  •  ^  higher  order  terms 


(2) 


n 
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The  approach  used  in  this  routine  is  a  first-order  approach 
and  neglects  the  higher  order  terms.  The  values  of  the  inde¬ 
pendent  variables,  x^,  X2,  ...  at  the  point  where  the 
function  is  to  be  evaluated  are  labeled  XV(1) ,  XV(2) ,  ...  XV(n) 
in  the  subroutine. 

The  subroutine  is  divided  into  three  main  sections  as 
shown  in  Figure  3.  The  first  step  in  applying  Eq  (2)  is  to 
search  the  data  to  find  the  indices  of  the  dependent  and  inde¬ 
pendent  variables  at  a  nearby  known  data  point  to  be  used  for 

(x,' ,  xl,  ...  x').  The  next  section  of  the  routine  checks  for 
1  2  n 

missing  values  (zeroes)  in  the  dependent  data  locations  speci¬ 
fied  above  and  increments  or  decrements  the  indices  to  find 
the  nearest,  non-zero,  dependent  variables  above  and  below  the 
test  point. 

In  the  final  section  of  the  routine  the  partial  differ¬ 
entials  are  approximated  by  first  difference  expressions  which 
are  evaluated  in  each  dimension  and  then  summed  to  find  the 
desired  output.  In  the  remainder  of  this  chapter,  these 
three  sections  are  discussed  in  more  detail. 


INTER? 


SEARCH 
FOR  INDICES 


CO^ffENSATION 
FOR  MISS’ING 
DATA 


DIFFERENCE 

EXPRESSICNS 


C 


RETURN 


Fig  3.  Flowchart  of  Main  Sections  of  INTER? 
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Finding  Index 

The  first  section  of  the  routine  contains  two  main 
loops,  each  with  an  internal  loop  as  seen  in  Figure  4.  The 
first  main  loop  searches  each  dimension,  or  string  of  indepen¬ 
dent  variable  values  (VAR) ,  to  find  the  index  of  the  lowest 
value  which  is  greater  than  or  equal  to  XV(I).  When  this  value 
is  found,  it  is  stored  in  LWK(I) .  If  XV(I)  is  less  than  the 
lowest  value  or  greater  than  the  highest  value  in  the  string, 
an  error  code  is  formed  and  printed,  and  control  returns  to 
the  main  program  (note  paths  A  and  B  in  Fig  4) .  The  second 
main  loop  of  this  section  uses  the  computed  values  in  LWK  and 
the  values  in  the  length  vector,  LV,  to  compute  INDEX.  INDEX 
corresponds  to  the  second  index  of  WVE  (the  known  dependent 
variable  array)  at  a  point  near  the  desired  value  specified 
in  XV. 

The  error  indicator,  lERR,  is  initialized  to  zero  at 
the  beginning  of  the  routine  and  is  revalued  when  necessary 
by  the  output  error  code  as  shown  in  Figure  5.  The  most  signi¬ 
ficant  digit  of  lERR  comes  from  I,  the  loop  index,  and  indicates 
the  dimension  being  searched  when  the  out-of-range  condition 
occurred.  The  least  significant  digit  is  left  zero  in  the  case 
of  a  below  range  condition  (path  A) ,  and  is  changed  to  nine 
in  the  case  of  a  value  of  XV  beyond  range  (path  B) . 

Compensation  for 
Missing  Data 

Section  2  of  the  routine  is  actually  within  two  do 
loops  of  section  3  (see  Figs  6  and  7) .  The  discussion  of  both 
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Fig  4.  Flowchart  of  Section  1,  Finding  Index 
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Fig  S.  Flowchart  o£  Output  Error  Code 


sections  will  be  clearer  i£  section  2  is  considered  £irst. 

It  must  be  recalled  that  TMP2,  ITl,  and  IT2  are  already 
assigned  in  section  3  and  that  section  2  is  included  within 
loops  with  indices  K  and  L. 

The  source  o£  the  data  used  as  inputs  to  this  routine 
are  usually  the  result  o£  point  designs  developed  by  the  Aero 
Propulsion  Laboratory  or  by  civilian  contractors.  Examining 
existing  data  reveals  that  in  many  cases  data  £or  some  o£  the 
points  within  the  array  are  not  available.  In  this  routine 
the  missing  values  o£  WVE  are  £illed  with  zeroes.  Since  these 
values  normally  represent  values  of  weight,  volume  or  e££ici- 
ency  o£  a  system  component,  a  value  o£  zero  should  not  other¬ 
wise  appear.  The  £inal  section  o£  the  program  will  use  the 
values  o£  the  dependent  variables  at  points  above  and  below 
the  test  point  in  each  dimension.  This  requires  that  many 
points  be  checked  £or  zero. 
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LHI=LHI+1 

INDHI*INDHI 

+IT1 


I=L 


See  Figure  7  for  C  and  D 


Fig  6.  Flowchart  of  Section  2, 
Compensation  for  Missing  Data 


Fig  7.  Flowchart  of  Section  3,  Difference  Expressions 
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In  the  first  part  of  this  section,  the  index  of  the 
point  above  the  test  point,  INDHI ,  is  initialized  as  index 


I 

(see  Fig  6).  The  stored  value  of  WVE  (K, INDHI)  is  then  tested 
for  zero.  If  a  zero  is  found,  INDHI  is  incremented  and  the 
check  is  made  again.  If  the  check  fails  after  the  highest 
value  in  the  string  has  been  checked,  an  error  output  is  gen¬ 
erated  (path  b)  as  in  section  1,  indicating  that  the  point  is 
beyond  range  and  in  which  dimension.  Note  that  the  index  of 
the  independent  variable,  originally  stored  in  LWK,  is  stored 
in  LHI  and  is  incremented  along  with  INDHI.  This  value  is 
required  in  the  final  section  of  the  routine. 

The  remainder  of  this  section  repeats  the  process  just 
described  for  the  known  point  below  the  test  point.  In  this 
case,  the  dependent  variable  index  is  stored  in  INDLO  and  the 
dependent  variable  index  in  LLO.  Again,  a  route  to  the  error 
code  is  provided  (path  A)  when  a  non- zero  value  is  not  found 
in  the  string. 

Difference  Expressions 

The  objective  of  the  final  section  is  to  determine  the 
estimated  values  of  the  dependent  variables  at  the  test  point. 
This  is  done  by  approximating  Eq  (2),  using  first  difference 
equations  to  replace  the  partial  differential.  TMPl  divided 
by  VAR(IT2+LHI)  -  VAR(IT2+LL0)  replaces  the  partial  differen¬ 
tial  and  (x^  =  x^)  is  replaced  by  XV(L)  -  VAR(IT2+LHI) .  Com¬ 
pare  Figure  7  and  Eq  (2).  The  quantity  WVE(K, INDHI)  in  the 
^  final  statement  corresponds  to  F(x^,  x^,  ...  x^)  in  Eq  (2). 
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The  results  of  these  approximations  are  discussed  in 
the  following  chapter.  It  is  shown  there  that  the  dependent 
variable  values  are  output  unchanged  when  the  test  point  coin¬ 
cides  with  a  known  data  point,  and  that  the  output  corresponds 
to  a  straight-line  approximation  for  points  between  known  data 
points. 
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V.  Evaluation  of  SUBROUTINE  INTERP 

The  objective  of  the  evaluation  was  to  choose  a  repre¬ 
sentative  data  base  and  develop  the  program  to  call  INTERP. 

The  output  of  this  program  is  then  studied  to  ensure  that  the 
subroutine  is  producing  reliable  data.  The  data  set  selected 
contains  the  results  of  a  conventional  generator  design  pro¬ 
gram. 

This  data  set  is  well  suited  for  a  test  of  the  routine 
with  three  independent  variables  and  three  dependent  variables. 
The  independent  variables  are  RPM,  voltage,  and  power.  The 
dependent  variables  are  specific  weight,  efficiency,  and 
volume.  This  data  set  has  the  added  advantage  that  an  algorithm 
has  already  been  developed  to  approximate  the  variation  of  spe¬ 
cific  weight  as  a  function  of  the  three  independent  variables. 
This  algorithm  is  the  one  presented  earlier  as  a  representative 
case  of  three  algorithms  (Eq  1,  page  7).  This  will  provide 
the  opportunity  to  compare  the  subroutine  output  with  the 
algorithm  results. 

Program  CGEN 

The  program  developed  to  manipulate  the  conventional 
generator  data  is  named  CGEN.  A  complete  listing  of  the  pro¬ 
gram,  along  with  input  and  output  data,  is  located  in  Appendix 
B.  As  seen  in  the  flowchart  in  Figure  8,  the  program  begins 
by  printing  the  labels  for  the  output  data.  Both  the  dependent 
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Fig  8.  Flowchart  of  Program  CGEN 


and  the  independent  variable  values  at  the  test  points  are 
included  in  the  printed  output.  The  dimensioning  information 
and  values  of  the  independent  and  dependent  variables  at  the 
known  points  are  then  read  from  the  data  file.  As  anticipated, 
there  were  numerous  cases  where  the  dependent  variable  values 
were  not  available.  Zero  fill  was  used  at  these  positions, 
as  previously  discussed. 

The  remainder  of  the  program  reads  a  test  point,  XV, 
and  calls  INTERP.  If  no  error  is  encountered,  the  output  is 
printed.  The  program  then  computes  the  specific  weight  using 
the  algorithm  presented  earlier  (Eq  1)  and  prints  that  value 
for  comparison.  If  an  error  is  encountered,  the  message  will 
be  printed  by  the  routine  and  the  program  goes  on  to  the  next 
test  point.  The  program  stops  when  it  has  used  all  the  test 
points  in  the  input  data  file. 

Results 

Many  test  points  were  run  through  CGEN  to  ensure  that 
the  routine  was  working  properly.  Several  errors  in  the 
routine  were  discovered  and  corrected  before  it  reached  its 
present  form.  Test  points  were  run  at  known  points  to  ensure 
the  routine  agreed  there.  Points  above  and  below  the  range 
were  used  to  test  the  error  code  function.  Points  were  also 
tried  at  and  near  locations  with  zero  fill  to  ensure  section 
2  of  the  routine  performs  correctly. 

With  all  corrections  made  and  all  functions  of  the 
program  and  routine  working  correctly,  the  test  points  included 


VOLTS  (L-N)  »10‘ 

Fig  9.  Comparison  of  Algorithm  and  Subroutine 
Changing  Voltage,  Volts  .65 


in  Appendix  B  were  used  to  produce  a  set  of  data  for  closer 
study.  In  one  group  of  test  points,  the  voltage  was  varied 
while  holding  the  RPM  at  18,000  and  the  power  level  at  8.83 
megawatts.  The  results  are  plotted  in  Figure  9.  As  can  be 
seen  in  the  figure,  the  output  from  INTERP  agree  with  the 
given  data  at  the  known  points  and  form  a  straight-line 
approximation  between  known  points.  The  value  from  the  algor¬ 
ithm  is  close  to  the  given  data  at  the  lowest  voltage,  an 
error  of  6.6  percent.  However,  at  the  highest  voltage  the 
error  increases  to  40  percent. 
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There  are  two  additional  groups  of  data  presented  in 
Appendix  B,  In  one  group  the  voltage  is  held  constant  at 
577  volts  and  the  power  at  4.91  megawatts.  The  RPM  is  varied 
from  8,000  to  20,000.  The  subroutine  output  again  agreed  with 
the  known  data  and  a  straight-line  approximation  between  points. 
The  algorithm  data  agreed  much  better  in  this  case  with  the 
error  varying  between  1.2  and  12  percent. 

In  the  other  portion  of  the  data  set,  the  RPM  was  held 
at  18,000  and  the  voltage  at  2885  volts,  while  the  power  level 
was  varied.  The  subroutine  output  again  agreed  well  with  the 
data  as  expected;  however,  the  error  in  the  algorithm  output 
was  nearly  constant  at  40  percent. 

Plots  of  all  three  sections  of  the  data  are  included 
in  Appendix  B.  Any  number  of  data  sets  and  plots  could  have 
been  listed  and  plotted,  but  these  three  were  chosen  to  show 
the  range  of  possibilities  and  advantages  of  the  subroutine 
approach  over  the  algorithm  approach  used  previously. 
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VI .  Conclusions  and  Recommendations 

The  original  objective  of  this  study  was  the  develop¬ 
ment  of  a  computer  program  to  produce  the  weight  and  volume 
algorithms  required  for  the  computer  aided  design  of  high 
power  systems  (Fig  1,  page  2).  These  algorithms  were  to  be 
based  on  the  data  contained  in  point  designs  being  produced 
by  contractors.  A  careful  comparison  of  the  advantages  and 
disadvantages  of  an  algorithm  program  versus  a  subroutine  to 
assess  and  interpolate  the  data  from  these  designs  was  made. 

The  result  of  that  comparison  was  the  decision  to  develop  an 
interpolating  subroutine. 

Conclusions 

INTER?  is  the  subroutine  that  was  developed.  It  has 
the  flexibility  to  interpolate  a  data  array  with  two  or  more 
independent  variables  and  will  output  the  interpolated  values 
of  any  number  of  dependent  variables.  Missing  values  in  the 
stored  data  are  compensated  for  by  using  zero  fill  at  these 
data  locations.  Output  values  at  known  points  agree  exactly 
with  the  stored  values,  and  the  values  between  points  are  based 
on  a  straight-line  approximation  in  each  dimension. 

The  results  of  the  evaluation  made  show  that  the  rou¬ 
tine  meets  all  the  objectives.  The  flexibility  of  the  routine 
exceeds  the  original  requirement  to  handle  five  independent 
and  three  dependent  variables.  Accuracy  is  as  expected,  and 
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an  error  code  informs  the  user  when  a  desired  test  point  is 
outside  of  the  data  range.  A  comparison  with  a  previously 
developed  algorithm  shows  that  INTERP  can  greatly  reduce 
errors  that  sometimes  exceeded  40  percent  using  the  algorithm 
approach. 

Recommendations 

One  of  the  next  steps  in  this  computer-aided  design 
program  will  be  the  development  of  an  executive  program  to 
implement  the  design  tradeoff  study  in  Figure  1.  This  program 
will  iterate  through  large  numbers  of  possible  component  con¬ 
figurations  to  find  designs  which  optimize  weight,  volume, 
and  efficiency. 

Several  of  the  component  design  programs  have  already 
been  developed  and  a  major  task  will  be  the  interfacing  of 
these  programs,  available  data  sets,  and  the  subroutine 
INTERP.  Another  project  that  may  prove  beneficial  is  the 
extension  of  INTERP  to  give  it  the  capability  to  also  extra¬ 
polate  beyond  the  given  data  range.  It  is  not  yet  clear 
whether  there  will  be  a  requirement  for  this  capability. 
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User  Guide  to  INTERP 

This  routine  is  used  to  interpolate  between  known 
data  points  within  a  data  array  containing  two  or  more  inde¬ 
pendent  variables  and  any  number  of  dependent  variables.  The 
user  must  supply  the  routine  with  known  values  of  both  the 
independent  variables  (VAR)  and  the  dependent  variables  (WVE) 
at  some  known  points.  Values  for  the  independent  variables 
are  then  input  (the  test  point)  and  the  outputs  of  the  routine 
are  the  dependent  variable  values  associated  with  the  test 
point. 

If  a  value  of  the  test  point  is  outside  the  range  of 
the  known  values,  an  error  message  will  be  output,  indicating 
which  dimension  (or  which  of  the  independent  variables)  was 
out  of  range,  and  whether  the  test  value  was  below  or  beyond 
the  data  range.  This  error  code  (lERR)  is  further  discussed 
in  the  argument  description  below. 

To  call  INTERP,  the  following  statement  is  used: 

CALL  INTERP  (NDIM,  LV,  LVT,  LDI ,  LDJ,  VAR,  WVE,  XV, 

LWK,  OUT,  I ERR) 

Arguments  must  be  dimensional  and  valued  as  described  below: 

NDIM  --  Number  of  dimensions  or  independent  variables 
(integer,  two  or  higher) 

LV  --  Length  vector  containing  number  of  known  data 
points  in  each  dimension  (integer  vector, 
length  equals  NDIM) 
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LVT 


LDI 

LDJ 

VAR 

WVE 

XV 

LWK 

OUT 

lERR 


Sum  of  the  elements  in  the  vector  LV  (integer) 
Number  of  dependent  variables  (any  positive 
integer) 

Product  of  the  elements  in  the  vector  LV 
(integer) 

Values  of  the  independent  variables  at  the 
known  points  (vector,  length  equals  LVT) 

Values  of  the  dependent  variables  at  the  known 
points  (array,  LDI  by  LDJ) ,  use  zero  fill  for 
array  points  with  no  data 

Test  point,  coordinates  where  the  dependent 
variables  are  to  be  estimated  (vector,  length 
equals  NDIM) 

Working  space  (vector,  length  equals  NDIM) 
Output,  containing  computed  values  of  the 
dependent  variables  (vector,  length  equals  LDI) 
Error  message  (output) ,  zero  indicates  no  error 
otherwise,  first  digit  indicates  dimension, 
second  digit  indicates  out  of  range;  zero  - 
below  data  range,  nine  -  above  data  range. 


Sample  Program 

Program  SAMPLE  is  listed  below  as  an  example  to  the 
user.  The  known  data  is  shown  in  Table  I.  Notice  there  are 
two  independent  variables  (VARl,  VAR2)  and  one  dependent  vari¬ 
able  (WVE).  If  there  was  a  second  dependent  variable,  its 
values  would  have  been  stored  in  WVE  (1,7)  through  WVE  (1,12). 


'! 


The  test  points  used  and  the  resulting  output  are  shown 
in  Table  I  for  comparison  with  Table  I.  If  one  of  the  dependent 
variable  values  listed  in  Table  I  had  not  been  available,  it 
could  have  been  replaced  by  a  zero.  INTER?  would  then  have  ig¬ 
nored  the  point  and  made  the  interpolation  based  on  the  remain¬ 
ing  points. 

TABLE  I 


Input  Data  for 

Program 

SAMPLE 

^^ARl 

VAR2V 

S.O 

10.0 

20.0 

20.0 

1.0 

2,0 

3.5 

50.0 

3.0 

4.0 

7.6 

WVE 

TABLE  II 


Test  Points  and  Resulting  Output  (SAMPLE) 


VARl 

VAR  2 

OUT 

6 

25 

1.53 

8 

.  30 

2.27 

12 

40 

3.35 

16 

45 

5,47 

18 

55 

* 

*Error  (IERR=29),  VAR2  Beyond  Range 
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PPOGPRM  xhMPLE  <  I  HPlJT .  OUTPUT') 

D I  MENS  I  ON  I.IVE  •:  1 .  t. )  »  VRP  <:  5>  •  XV  <  £ )  .  LMH:  »  OUT  C 1 ")  » LV  < 

C 

C  II I  MENS  I OU  IHFOPMRTION 
C 

MriIM=£ 

LV<n=3 
LV<£>=d 
LVT=5 
LDI  =  1 
LriJ=6 
C 

C  KNOl.lU  lUriEPEHIiENT  VhPIRBLE  VRLUES 

c 


VRR 

1 

=5 

VRR 

=  1 

0 

VRR 

■  3> 

0 

VRR 

:4) 

=2 

0 

VRR 

5) 

0 

c 

C  KNOl.lM  IIEPENDEHT  VRPIRELE  VRLUES 
C 

1..1VE  -'It.  r)=i 
lilVE  '  1  !•  c! )  — £ 

I.IVEa»3)=3.5 
I.IVE  <;  1 »  4 )  =3 
1,1  VE  <  1 1  5  :■  =4 
I.IVE  Cl*  6)  =7. 6 
C 

C  PERU  TEST  POINT 
15  P  E  R  Ii  ♦ '  •  X V  .  !  ':•  »  I  =  1 » NH I M') 

I F  ■:  EOF  .:  5L I NPUT )  .  NE .  0 .  >  STOP 
CRLL  INTEPP  CNIiir-U  LV>  LV'T?  LDI ..  LUJ.  VRP*  I.IVE » 
1  XV<Ll.il:..DUT?  IERP> 

PRINT**  ''Ol.lT  =  ‘'*OUTa') 

GO  TO  15 
END 


Input  Data  Output  Data 


6  £5 
8  3  0 

12  40 
16  45 

13  55 
♦EDR 
♦EOF 


0:jT=£.  36666 
□UT=:;;.  35333 
0*,lT“5. 4766'6 
IE.RR=£9 
□UT=5.47666i 
STOP 
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Listing  of  INTER? 


SUBF'nUTINE  INTEPP  f  NIiIM!.  LVj  LVT»  LDI  >  LDJ>  VPP»  WVE>  XV, 
1  L.:.i!;  ,a!jT,  lEPP.) 

D I MEMS I  MVE  ••  LD I ,  L’j  J;-  ,  V:=lR  -rLVT  :-  ,  XV  >: ND I  M>  , 

1  LMi<  •  NniMJ  ,  nUT  (LDI.;'  ,  LV  CNDIM.) 

’ECTinM  1,  FINDir^G  INDEX 


ICi 


80 


30 

31 

35 


IEPP=0 
MIN=1 
MPX^LV  ( 1  :• 

FIND  LVir?  !□  FEPM  INDEX 
Dd  80  I=1,NDIN 

CHECK  FOP  VPLLIE  BECaM  PHNGE!  IF  SO,  GC  Tn  EPPDP  CDDE 

IF(XVa>  .LT.  VflP  (:MIN>>  go  TD  30 
LST=LV<I>  fl 

L.nnP  TO  FIND  THE  MININUN  STDPED  VHLUE,  WHICH  IS  GPEHTEP 
THHN  CP  EQiJHL  TO  THE  TEST  POINT 

DO  5  K=1,LST 
J=K-1 

I F  -rXV  I  >  .  LE .  VHP  CM  I N  J  )  >  GO  TO  1  0 

CHECK  FOP  VHI-iJE  PEOVE  PRNGE;  IF  SO,  GO  TO  ERPOP  CODE 

IFCXVCD  .GT.VPPCMPX^  >  GO  TO  35 
5  CONTINUE 

STOPE  THE  INDEX  OF  THE  INDEPENDENT  VHPIRELE  FOUND  IN  LWKCI) 
!_:.IK  (1:5=1  t-J 

SET  MPX  PND  MIN  FOP  NEXT  DIMENSION 

MRX=MPX<-LV  Cl  ♦■i;5 
MIN=MIN<-LVCI> 

GO  TO  39 

OUTPUT  EPPOP  CODE 
IEPP=I»10 

PPINT», "lEPP*", lEPR 
PETUPN 
IEPP=I»10«-9 
GO  TO  31 
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o  o  a  n  n  a  n  o  o  a  n  o  o  n  n 


C  riNP  INDEX  GF  MVE  CGEFEGFONDINF  TG  LMK 

39  INDEX=0 

DG  50  N=E'>NDIM 

ITEMP=1 

NDE=N-1 

DG  40  M=1.NDE 

4  0  I TEMF  - 1  TEMP  V  ''M.) 

■50  rNriEX=  •:LMK  <N>  -1  >  ♦ITEMP^-INDEX 

INDE>'=INDEX4-LI.I}<:  a> 

SECTIGN  3.  LGGPS  TG  FGPM  GUT  ■:K> 

DG  70  K  =  l!.|_DI 
TMP£=0 
ITE=0 
IT1  =  1 

DG  SO  L=1jNDIM 

SECTIGN  cN  CGMPENS9TIGN  FGP  MISSING  DPTft 

INITIALIZE  INDUI  S-CHECK  FGP  ZEPG  AT  MVE  •:!<  >  INDHI> 

IND!-II  =  INDEX 
LHI=LMI::  CL) 

MMAX=LVCL>  -L1..IK<;L> 

DG  55  M=1.MMAX 

I F  c  !.l VE  c  l<  .  I NDH I  >  .  NE .  0>  GG  TG  5S 
IF  ZEPG»  INCPEMENT  I NDH I  AND  P  'ECK 
LHI=LHIfl 

55  INDHI  =  INDHH-IT1 

IF  END  GF  DATA  STPING  IS  PEACHED^  SEND  EPPGP  MESSAGE 


V# 


I=L 

GG  TG  35 

INITIALIZE  INDLG  e  CHECL  FGP  ZEPG  AT  MVE  ck .  INDLG;) 

5S  INDLG=INDEX-IT1 

LLG=LM!-  CL  ) -1 

NMAX=LLa 

DG  57  N=1,NMAX 

IF  c  (MVE  (L.  INDLG;)  .NE.  0.)  .'GP.  c;xv<L.)  .EC.  VAP  c  ITE^-LHI  ) 
GG  TG 

IF  ZEPG  DECPEMENT  INDLG 
LLG=LLa-l 

57  INDLG=INDLa-ITl 
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TABLE  III 


Input  Data  for  Program  CGEN 


*No  Data  Available 


Independent  Variables 

RPM 

(*  1000) 

Voltage 

(volts) 

Power 

(^nv) 

2308 

2308 

2308 

2308 

2308 

2308 

2308 

2885 

2885 

2885 

2885 

2885 

2885 

2885 


4.91 

4.91 

4.91 

4.91 

4.91 

4.91 

4.91 

4.91 

4.91 

k  r\  1 


4.91 

4.91 

6.87 

6.87 

6.87 

6.87 

6.87 

6.87 

6.87 


Dependent  Variables 


Sp.  Wt.  Efficiency 
(Ib/kW)  (%)  *1000 


17.3 
10.6 
7.6 
6.0 
5.1 
.6 
.2 


94.7 

95.4 

95.6 

95.7 
95.7 

95.5 

95.3 

94.6 

95.4 

95.7 

95.7 

95.5 

95.5 

95.3 

95.4 

95.8 
95.8 

95.6 

95.3 
95.0 

94.3 


17.6 

10.7 
7.8 
6.2 
5.1 

4.7 
4.5 

17.7 
11.3 

8.4 

6.8 
6.0 
5.3 
4.9 


Independent  Variables  j 

RPM 

C*  1000) 

Voltage 

(volts) 

Power 

(MW) 

8 

1154 

6.87 

10 

1154 

6.87 

12 

1154 

6.87 

14 

1154 

6.87 

16 

1154 

6.87 

18 

1154 

6.87 

20 

1154 

6.87 

8 

1731 

6.87 

10 

1731 

6.87 

12 

1731 

6.87 

14 

1731 

6.87 

16 

1731 

6.87 

18 

1731 

6.87 

20 

1731 

6.87 

8 

2308 

6.87 

10 

2308 

6.87 

12 

2308 

6.87 

14 

2308 

6.87 

16 

2308 

6.87 

18 

2308 

6.87 

20 

2308 

6.87 

Dependent  Variables 


Sp.  Wt.  EfficiencJ 
(Ib/kW)  I  I 


95.5 

95.8 
96.0 
96.0 

95.6 

95.5 

* 

95.6 

95.7 
96.0 

95.8 

95.7 

95.5 

* 

* 

95.8 
96.0 

95.8 
95.7 

95.6 
* 


11.5 

8.6 

6.9 

6.2 

5.4 

* 

18.2 

11.8 

8.7 

7.2 

6.3 

5.7 

* 

* 

12.1 

8.9 


*No  Data  Available 


Independent  Variables 


Dependent  Variables 


RPM 

(*  1000) 

Voltage 

(volts) 

Power 
(Mv; ) 

Sp .  Wt . 
(Ib/kW) 

E££icienc> 

(%) 

8 

2885 

6.87 

.  215 

95.4 

18.8 

10 

2885 

6.87 

.186 

95.8 

12.2 

12 

2885 

6.87 

.173 

95.8 

9.3 

14 

2885 

1 

6.87 

.154 

•  96.0 

7.6 

16 

2885 

6.87 

S  .151 

95.6 

6.7 

18 

2885 

6.87 

.142 

95.5 

6.0 

20 

2885 

6.87 

* 

1 

* 

* 

8 

1 

1154 

8.83 

.188 

95.7 

19.1 

10 

1154 

8.83 

.168 

96.0 

12.8 

12 

1154 

8.83 

.153 

95.9 

9.9 

14 

1154 

8.83 

.136 

95.9 

8.1 

16 

1154 

8.83 

.138 

95.4 

7.5 

18 

1154 

8.83 

.131 

95.1 

6.8 

20 

1154 

8.83 

* 

* 

* 

8 

577 

8.83 

* 

* 

* 

10 

577 

8.83 

* 

* 

ft 

12 

577 

§.83 

* 

* 

* 

14 

577 

1 

8.83 

* 

* 

* 

16 

577 

8.83  ' 

* 

* 

* 

18 

577 

8.83 

* 

* 

* 

20 

577 

8.83 

* 

* 

* 

*No  Data  Available 


o  o  o  n  o  o  o  o  o  o  o  o  o  n  o  o  o  o  non  o  o  o 


Listing  of  Program  CGEN 

PPnGPHM  r  I3EH  ■  I NPUT  .  OUTPUT  ' 

D I  f '1 E  r  1  i  1 0 1  ^  1 1 1 V  E  '  !• '  1 1  •  5  >  »  V  H  P  1 15  ■'  «  X  V  '•  3  ■'  •  L I.'  I K  '  3  >  <0  U  T  ••  3  '  •  L  V  '•  3  ■' 

C 

C  PRINT  COLUMN  HEftlUNGS  FDR  DUTRIJT  POTfl 
C 

£  PRINTS.  "RPN  ♦K  V  L-N  PWP  CMW;'  SPMT  LE/Kl,.j  EFF  V 

IIN" 

4  FDRNPT  (£  (F5.  0*5X)  .4'rF6.  3!.6X)> 

PEflli  INFDPMPYIDN  NEEDED  TO  DIMENSION  ARRAYS 

READ^-  NDIN.LDI  ?LDJ!.LVT 
READ^'  'LV' :l> • I  =  1»NDIM) 

READ  THE  INDEPENDENT  VARIABLE  VALUES  ‘ 

PE AD^  •  '•  VAR  I  I  j  »  1  =  1*  LVT.> 

READ  THE  DEPENDENT  VARIABLE  VALUES 

DO  10  J=l!.LDI 

1  0  READ^»  tl'lVE  ':..Uk  )  »K=1»  105> 

READ  THE  TEST  POINT 

1 5  READ^ »  ■:  XV  I  >  »  I  =  1 » ND IM> 

IF '  EOF 'iS'LINPUT':' .NE.  0.  :>  STOP  "END  OF  RIJN" 

CALL  I NTERP  •:  ND  I M .  LV » LVT .  LD I  ■«  LD  J .  VAR » I.IVE « 

1  YV.LUK -OUT, IEPR> 

IF  THERE  I'.lAS  AN  ERROR?  GO  TO  THE  NEXT  POINT 
IF'IEPR.NE. 0>  GO  TO  15 

PRINT  THE  INDEPENDENT  AND  DEPENDENT  VALUES  AT  THE  TEST  POINT 

PRINT  4.  (  cXV':!:)  <  I  =  l?3>  ?  ■rOUT':  J>  » .J=l?3>> 

COMPUTE  THE  SPEC.  I-IT.  USING  THE  ALGORITHM?  FDR  COMPARISON 

S  l.i  1  = .  1 5  7  ♦  <'  1 .  £  9  - .  £  9  ♦  <’  X  V  3  :•  •  ••  5  >  ♦♦.449'’  ♦  <  - ,  09+1 .  0  9  ♦  <  X  V  >:  1  >  / 1 4  > 

►  ♦  9£  05  '  'j  )  ♦  ':.  9597+.  1433+  '^XV  •:£>  ♦SORT  (3.  >  -- 1  00 0?  > 

PRINT  THE  SPEC.  I.IT .  COMPUTED 

PR  I  NT  ♦  ?  "  SP  I.IT  ■:  FROM  ALG  >  =  "  ?  SW 

GO  TO  15 

END 


J 
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Output  from  CGEN 

1  05  •  15 
7  5  3 

10-  <  14.  16.  18.80. 5.77.  1154.  1731.8308. 


4.91.6.87.8.8  3 


■  c'c'4 

.  194 

.160  .140  .134  .184 

.  188 

.  831 

.185  .164  .145  .138 

.  188 

.  138 

.886  .195  .168  .157 

.  146 

.  133 

.000  .^37  .195  .176 

.  16  0 

.  149 

.137  .134  .845  .197 

.  170 

.  1 6  0 

. 149  . 147  . 148  . 196 

.  178 

.  154 

.139  .137  .181  .180 

.  197 

.  175 

.154  .140  .134  .185 

.  0  00 

.  801 

.188  .156  .149  .138 

.  1 3  0 

.  0  0  0 

.000  .177  .160  .154 

.  1 43 

.  135 

.000  .815  .186  .173 

.  154 

.  151 

.148  . 000  . 000  . 000 

.  0  0  0 

.  00  0 

.000  .000  .000  .188 

•  1 6<y 

.  153 

.  136 

.138  .131  .000  .185 

•  1 66 

.  158 

.148  .141  .188  .000 

.  1  S‘? 

.  171 

.156  .151  .145  .141 

.  000 

.  1 95 

.173  .158  .154  .147 

.  1 37 

.  0  0  0 

94.9 

95.  3 

95.5  95.8  95.4  95.4 

94.  9 

94 . 8 

95.5  95.7  95.8  95.7 

95 . 6 

95.  0 

94.9  95.4  95.8  95.6 

95.6 

95 . 6 

0  0 . 0  94 . 7  95 .4  95.6 

95 . 7 

95 . 7 

95.5  95.3  94.6  95.4 

95.7 

95. 7 

95; 7  95.5  95.3  95.4 

95.8 

95.8 

95 . 6  95 . 3  95 . 0  94 . 3 

95.5 

95.8 

96. 0  96. 0  95.6  95.5 

0  0 .  0 

95.6 

95.7  96.0  95.8  95.7 

95 . 5 

0  0.  0 

00.0  95.8  96.0  95.8 

95 . 7 

95.  6 

0  0 .  0  95 . 4  95 . 8  95 . 8 

96 .  0 

95.6 

95.5  0  0.  0  00.  0  00.0 

0  0 .  0 

0  0 .  0 

0  0 .  0  0  0 .  0  00.  0 

95.7 

96.  0 

95.9  95.9  95.4  95.1 

00.  0 

35.  S 

96.  1 

36  •  1 

‘^5  .  '^5 . 6.  .  5  0  0 ,  0 

35 . 3 

96.  0 

96.1  95.8  95.6  95.8 

00.  0 

95.8 

96.0  96.1  95.8  95.7 

95 . 5 

00.  0 

16.7 

10.  1 

7.0  5.5  4.8  4.1  3.8 

16*3 

10.8 

7.8  5.6  4.9  4.8  4.0 

17.  0 

10.5 

7.4  5.9  5.  0  4.4  0 .  0 

17.3 

1  0 .  6 

i‘  *6  t*  •  •  1  *4  ■  6  ■  ib 

17.6 

10.7 

7.8  6.8  5.1  4.7  4.5 

17.7 

11.3 

8 .4  6 . 8  .0  5.3  4.9 

18.  0 

11.5 

8 .6  6.9  6.8  5.4  0 .  0 

18.8 

11.8 

•r*  m  I’  i’  •  c'  6  •  •  0 

0.0  1 

8.1  8 

.9  7.4  6.4  5.9  0.0 

18.8  1 

8.8  9 

.3  7.6  6.7  6.  0  0.0 

0. 0  i; 

.  0  0 . 

0  0.  0  0.  0  0.  0  0.0 

19.  1 

18.8 

9.9  8.1  7 . 5  6 . 8  0 . 0 

19.5 

13.  0 

1 1.1  •  1  y  •  *4  r‘  •  i  '  6  •  S  0  ■ 

19.7 

13.3 

10.3  8.8  7.9  7.3  0. 

0  80.  1 

13.4 

10.5  9.0  8.1  7.3  0 

8  577  4.91 

9  577  4.91 


10 

577  4.91 

11 

577  4.91 

18 

577  4.91 

13 

577  4.91 

14 

577  4.91 

15 

577  4.91 

16 

577  4.91 

18  1154 

;3 . 9.3 

17 

577  4.91 

1 8  1 8  0  n 

8 . 83 

18 

577  4.91 

18  14  0  0 

8.83 

19 

577  4.91 

1 8  1 6  0  0 

8 .  S3 

80 

577  4.91 

'  18  1731 

8.83 

IS 

8885  4.91 

1 8  1 8  0  0 

8.83 

18 

8885  5. 

1 8  8  0  0  0 

y  ■  3  y 

IS 

8885  5.5 

18  8800 

3  •  3  3 

IS 

8885  6. 

1  y  c*  3  0  8 

8.83 

18 

8885  6.5 

18  8400 

8.83 

18 

8885  6.87 

1 8  86  0  0 

8 . 83 

18 

8885  7 . 

IS  8800 

8.83 

18 

8885  7.5 

18  8885 

8.83 

IS 

8885  8. 

♦EDR 

18 

8885  8.5 

♦EOF 

IS 

8885  8.83 

t 


I 

If 

I 

If 


't 


f 


\ 


Output  from 

CGEN 

PPM  ♦K 

V  !_-M 

PMP  CMf.l)  SPMT  LE-^-Kl.l 

ccc  i..'n> 

Cl  r  viJi_ 

KC'J  IN 

y . 

577. 

4.9iri 

.  ££4 

'94.  900 

18.700 

SP  MT 

'FPHM  :hL:5> 

=. ££85851 8583SS 

9. 

577. 

4.91  0 

.£•09 

95. 100 

13.400 

SP  :.iT 

CFPDM  RL:5:> 

= . £  0995  08545953 

1  0. 

577. 

4.910 

.  194 

95.300 

10.  100 

SP  MT 

<FPCM  :=iL:2;> 

=. 19808847994£ 

11. 

-•>1  1  ■ 

4.910 

.  177 

95 .400 

8.550 

SP  MT 

(FPDM  RL5:> 

~.  1  :z:4c'ir'5c’ M  U9'r*85 

1  c! . 

577. 

4.910 

.180 

95. 500 

7.  000 

SP  MT 

CFPDM  RL:3:) 

=. 174084££74S73 

IS-:. 

577 . 

4.91  0 

.  150 

•95.85  0 

8.  £50 

SP  MT 

CFPDM  RLG> 

=  .  18519l:54S080£ 

14. 

»■  f'  m 

4.910 

.  140 

95. 800 

5.500 

SP  MT 

CFPDM  RL:3:' 

=. 15734338£8831 

15. 

577. 

4.910 

.  137 

95.800 

5.  150 

SP  MT 

'•  F  F  n  M  JH  L  o  ■' 

=  .15  0354  0££59  09 

IS. 

577. 

4.91  0 

.  154 

95. 400 

4.800 

SP  MT 

cFPDM  RL:5:' 

=. 1440813001398 

17. 

577. 

4.91  0 

.  1  £9 

R5.400 

4.450 

SP  MT 

i:FPDM  RL:3> 

=  .  13:54134517£09 

18. 

577. 

4.910 

.  1£'4 

95.400 

4.  100 

SP  MT 

CFPDM  R1.8> 

=. 133£8 144479 08 

19. 

.'II. 

4.910 

.  1£3 

•95.  150 

3.  950 

SP  MT 

CFPDM  Rl_:3::' 

=. 1 £8553391 1983 

80. 

577 , 

4.910 

.  1££’ 

94. 900 

3.800 

SP  MT 

CFPDM  RL!3> 

=  .  1  £4£*3  05484  0 1 8 

18. 

OOOC 

4.91  0 

.  147 

95.500 

4.700 

SP  MT 

CFPDM  RL;3> 

= .  £  098  07£  0  04393 

18. 

£885. 

5.  000 

.  147 

95 .500 

4 .  I"’  8  0 

SP  1-1 T 

CFPDM  RDS::' 

=.£091315750777 

13. 

£885. 

5.500 

.  145 

•95.500 

5.  091 

SP  MT 

c  F  F  3  y  1_  o  -  ' 

=  .£08571  £"5583  01 

18. 

cyy5  • 

8.  000 

.  144 

'95 .  5  0  0 

5.4£3 

SP  MT 

■:FPDM  RL;3> 

=.£04138£99514S 

18. 

£885. 

S .  5  0  0 

.  143 

•95.500 

5.  755 

SP  MT 

CFPDM  RDG:- 

=.£018107499889 

18. 

£8:55. 

8.  870 

.  14£ 

•95.500 

8 .  0  0  0 

SP  MT 

CFPDM  RL:3.' 

=  .£•001538730145 

IS. 

£885. 

7.  OOm 

.  14£ 

•95.500 

8.  o:58 

SP  MT 

CFPDM  RL.S> 

=  .  1  •9-958 1^£74473 

IS. 

£885. 

7 .  5  0  0 

.  140 

'95 .  5  0  0 

8.418 

SP  MT 

'  FPDM  RL:’' 

=. 1974390187889 

IS. 

£885. 

8.  000 

.  139 

95.5  0  0 

8.749 

SP  MT 

<FPDM  RL:3) 

=. 19537358£3998 

IS. 

£:5:55. 

8.  500 

.  1 38 

95.500 

7.  081 

SP  MT 

CFPDM  RLG> 

=. 1933780834534 

IS. 

£885. 

8.830 

.  137 

95.500 

7.300 

SP  MT 

CFPDM  RLG> 

=. 19£ 0963 088485 
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RP“  * 

X  i-r 

OVF  (►'!■»)  SOMT  LB/<M 

£FF 

VHL  KC'l  T'J 

1ft. 

1154. 

8.  Sft  Cl 

.  131 

•95.100 

8.800 

:5ft  MT 

.:FftnM  !=iLft>=. 

Ift'ftft.ft 05 175885 

18. 

1  ft  Cl  0 . 

ft.  SftO 

.  131 

95. 13ft 

8 .  y  0  Cl 

:5ft  MT 

(FftHM  FlLS::'=. 

1410150215183 

IS. 

1 4  0  0 . 

ft.  83 0 

.  130 

95.271 

8.  .800 

5ft  MT 

(FftHM  ;=iL5'=. 

14707SO:Bft0'314 

1  ft. 

1  ft  0  0 . 

ft.  830 

.  1 29 

95.409 

8.800 

5P  MT 

■:FftnM  FiL5>=. 

1531411428845 

18. 

1 7.5 1 . 

ft.  830 

.  IftS 

95. 500 

8.800 

5ft  MT 

.:FftnM  ftLS>=. 

1571  lft'44733'38 

1ft. 

1  ft  Cl  0 . 

8.830 

.  130 

,  95.484 

8. 880 

:5P  MT 

■FftuM  ftLo;'=. 

1 5'3£'  04ft  03£:375 

IS. 

ft  0  Cl  0 . 

8.S30 

.  134 

95.380 

7.  033 

5ft  MT 

(FftHM  ftLft:'=. 

1 85£'87ft83:5 1  08 

1  ft'. 

ftftCiO. 

ft.  830 

.  1 39 

95.258 

7.206 

5P  MT 

(FftOM  :HLft>  =  . 

1713303243837 

1ft. 

ftftCift. 

8 . 83  0 

.  141 

95.200 

7.300 

5P  MT 

(FftaM  ftLG:'=. 

174804377093ft 

IS. 

ft  4  00. 

8.830 

.  140 

•95.  £48 

7. 

Sft  MT 

(FftOM  Fi!_8>=. 

1 773933849588 

18. 

ftftOO. 

8.830 

.  1 3'3 

95.  :35ft 

7.30. 

5P  MT 

cFfta.M  flL8:'=. 

1 834584455299 

18. 

£8  0  0 . 

8.  830 

.  1 38 

•95.456 

7 .  ft  L>  u 

5ft  MT 
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